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Abstract 

The random phase approximation (RPA) theory is used to derive the Luttinger parameter g for 
metallic carbon nanotubes. The results are consistent with the Tomonaga-Luttinger models. All 
metallic carbon nanotubes, regardless if they are armchair tubes, zigzag tubes, or chiral tubes, 
should have the same Luttinger parameter g. However, a (10,10) carbon peapod should have a 
smaller g value than a (10,10) carbon nanotube. Changing the Fermi level by applying a gate 
voltage has only a second order effect on the g value. RPA theory is a valid approach to calculate 
plasmon energy in carbon nanotube systems, regardless if the ground state is a Luttinger liquid or 
Fermi liquid. (This paper was published in PRB 66, 193405 (2002). However, Eqs. (6), (9), and 
(19) were misprinted there.) 
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In the past decade, the study of Luttinger liquids 1 has been fuelled by experimental ad- 
vances in one-dimensional (ID) systems such as semiconductor quantum wires, 2 and more 
recently, carbon nanotubes. 3 The non-Fermi liquid behavior of an interacting ID electron 
gas is predicted by Tomonaga-Luttinger models. 4 The assumed linear single-particle energy 
spectra render these models exactly solvable by bosonization method, yielding eigen modes 
as collective charge density modes (plasmons) and spin density modes. Single-particle exci- 
tations are forbidden in the Luttinger liquid. For wavevector q > 0, the single-particle exci- 
tation energy (for the non-interacting electron gas) in these models is linear in q, E s = qhvp, 
where vf is the Fermi velocity, while the collective charge density mode has energy 



/ IV 

Here V q is the eletron-electron interaction matrix element, and typically the value at q = 
is taken. The Luttinger parameter g is defined as the ratio of the two energies, 



This parameter is of key importance in the study of Luttinger liquids. It uniquely determines 
the exponents a of the power law behavior of the conductance with respect to temperature 
and voltage, and nonzero a values for carbon nanotubes are indications of a Luttinger 
liquid. 5 - 9 

In 1992, Li, Das Sarma, and Joynt 10 noted that the standard random phase approximation 
(RPA) theory applied to a ID semiconductor quantum wire produced results in quantitative 
agreement with the plasmon dispersion measured by Goni et al. 2 They also obtained an 
analytical solution for the plasmon energy in the case where only one subband is occupied 
by electrons. If the small q limit is taken, the analytical solution reduces to the result in 
Eq. (1). The agreement of the RPA theory with the exact result for Tomonaga-Luttinger 
models in the small q limit dispelled the myth that RPA was not suitable in ID. The need 
to take the small q limit in order to obtain agreement is due to the fact that in Ref. 10, a 
parabolic single-particle energy spectrum is used. Near the Fermi surface, the single-particle 
excitation energy for a non-interacting electron gas is 

+ g )'_^g = gte , + W (3, 
2m* 2m* H 2m* y ' 

where the Fermi wavevector kp is related to the Fermi velocity vp by vp = ftkpjm* . Only 

when q is small, is the last term in Eq. (3) negligible and the linear energy spectra in 



the Tomonaga-Luttinger models become good approximations. In real systems, the single- 
particle excitation energy is usually nonlinear. Hence the bosonization method and the RPA 
theory each has its advantages and shortcomings: The bosonization method is robust and 
exact for linear single-particle energy spectra, but is applicable only in the small q region 
where the dispersion is linear; The RPA theory is approximate in general but its validity is 
not limited to the small q limit. In RPA, the plasmon energy is determined by the zero of 
the dielectric function 

e(q,u) = l-V q P(q,u) = 0, (4) 
where P(q, uj) is the polarization of the non-interacting electron system, 

P(q, u>) = £ Pfa ^^E f kk+ \ kk h ■ (5) 
i L> i)k ^i,k+ q - &i,k - nu 

In the above equation, L is the length of the ID system, the factor of 2 comes from spin 
degeneracy, i is the subband index, E i k is the single-particle energy, and f it k is the Fermi 
distribution function. 

In this paper we derive the Luttinger parameter g for metallic carbon nanotubes. As is 
well known, the structure of carbon nanotubes is specified by an integer index pair (n, m) . 
Those with m = n are called armchair carbon nanotubes, those with m = are called zigzag 
carbon nanotubes, and others are called chiral carbon nanotubes. 11-14 All armchair carbon 
nanotubes are metallic, zigzag tubes with n being multiples of 3 are metallic, and chiral tubes 
with In + m being a multiple of 3 are metallic. Unlike in Ref. 10 where there is only one 
occupied subband and the single-particle energy dispersion is parabolic, carbon nanotubes 
have more (partially) occupied subbands and have sinusoidal band structures. 11-15 The band 
structure of carbon nanotubes can be easily derived from that of the graphite sheet, given 
by 11 " 13 

E graphlte (k x , ky) = ± 7 o[l + AcosC^yosi^f) + W(^)] 1 / 2 , (6) 

where a = 0.246 nm is the lattice constant of a graphite sheet, and 7 is the overlap integral 
between nearest-neighbor carbon atoms. By rolling up a graphite sheet along different 
directions, one obtains differently structured carbon nanotubes. 

The Brillouin zone of the graphite sheet is hexagonal, and k x and k y can be any value 
within the first Brillouin zone. Once the graphite sheet is rolled up into a nanotube, how- 
ever, allowed values for k x and k y are reduced due to quantization in the circumferential 
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direction. For the armchair carbon nanotube, k x is in the circumferential direction and can 
only take on quantized values k x = u2n / n^/3a, with v — 0, 1,2, ..,n. Most of the subbands 
in an armchair carbon nanotube are either fully below or fully above the Fermi energy 
Ep, with only two subbands crosssing the Fermi level (those obtained when v = n). 11-15 
For low energy (intrasubband) excitations, only these two subbands need to be considered. 
This is not an approximation, because other fully occupied or fully empty subbands only 
contribute to intersubband excitations, which have different angular momenta from intra- 
subband excitations. (In the carbon nanotube, the angular momentum is a good quantum 
number, therefore intersubband and intrasubband excitations are decoupled.) The energy 
dispersions of the two subbands crosssing the Fermi level are 

E hk = ± lQ [2cos(k y a/2) - 1]. (7) 

In Eq. (7) we take the convention that % =1 {% — 2) corresponds to the — (+) sign. The 
Fermi wavevector is determined by letting E iik — E F — 0, hence k F = 2ir/3a. Figure 1 
shows the dispersions of the two subbands. 

Near the K points where the two subbands cross the Fermi level and each other, the 
single-particle excitation energies are 

r-i r-i r-i , r + Q k F , . qd . ,k F d QCL v . . 

E s = E iikF+q - E iikF = ±27o[cos — - — a - cos —a\ = ^A>yQSin—sin(— + — ). (8) 

In the small q limit, to first order in q, Eq. (8) reduces to 

. k F a \/3 
E s = Tq7oasin{—) = TQ— 7oa. (9) 

If we compare this to the linear energy spectra of the Tomonaga-Luttinger models E s = qhvp, 
we see that in the case of armchair carbon nanotubes, we should equate 

hv F = 7 asm(^) = ^7oa. (10) 

Before we treat a real armchair carbon nanotube, let us first obtain the results for the 
case where there is only one subband crossing the Fermi level, i = 1 in Eq. (7). At zero 
temperature, the polarization is given by the integral 

PM = -[ (p p 1 Mh + r it tt^ ( n ) 
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where the range of integration is [— k F , k F ]. We obtain 

p 2 [hutan{k F a/A + qa/8)] 2 - [4 lo sin{qa / '4) + P^] 2 

, [hutan(k F a/4: - qa/8)] 2 - [4ry sin(qa/A) - PJ 2 

n [hutan(k F a/A - qa/8)} 2 - [4 7o sm(ga/4) + EJ\ 2 ' ' 1 J 

where P w = y / [47 sm(ga/4)] 2 — (ftu;) 2 . Note that E u is allowed to be imaginary. When q is 
small, tan(k F a/4:±qa/8) = tan(7c/6±qa/8) ~ l/\/3± ga/6. Doing Taylor series expansion 
to lowest order in Eq. (12) gives 

d/ \ 47 (ga) 2 

Pi (9, a;) = p-r 1 . (13) 

aiT\ / 3[^(huj) 2 - (7 ga) 2 ] 

To see how Eq. (13) is derived, note that when q is small the first log term in Eq.(12) 

is approximately (2/3v / 3)(^) 2 ga/[2(47 sm(ga/4)) 2 - (A/3)(huj) 2 - 8~f sin{qa/4)Ej[, while 

the second log term is approximately — (2/3A/3)(^a;) 2 ga/[2(47 sm(ga/4)) 2 — (A/3)(hu) 2 + 

8josin(qa/ '4)EJ\. Also the product of the denominators in these two expressions is exactly 

equal to (4/3)(huj) 2 [(4/3)(huj) 2 - (47 sm(ga/4)) 2 ]. From Eqs. (4) and (13), we can solve for 

the plasmon energy and obtain 

hw = q— 7004/ 1+ * ■ (14) 

We note that Eq. (14) is identical to Eq. (1) obtained from Tomonaga-Luttinger models if 
the equivalence in Eq. (10) is taken into account. 

Of course, in real armchair carbon nanotubes, we have to take into account both subbands 
in Eq. (7). For % = 2, the calculation of the polarization using Eq. (11) requires special 
care at the Brillouin zone boundaries. As we can see from Fig. 1, at zero temperature the 
i — 2 subband is occupied by electrons in two segments of the subband, [— n/a, —k F ] and 
[k F ,n/a\. However, if we let the integral of Eq. (11) to be over these two ranges, we get 
the unphysical result that the plasmon can exist only for strong and weak electron-electron 
interactions, but not for intermediate interaction strength. The cause of the problem is that 
in Eq. (11), we not only have subband energy E^, but also Ei^+q and Ei^-q, and the 
choice of the ranges [—n/a, —k F ] and [k F ,7i/a] did not ensure that k±q (q > 0) be confined 
within the Brillouin zone. The proper integration ranges for Eq. (11) are [—n/a,—k F ] 
and [k F ,n/a — q] for the first integrand, and [— n/a + q, —k F ] and [k F ,ir/a] for the second 
integrand. With the proper integration ranges, we find the two subbands to have the same 
polarization, P 2 (q,u) = P\{q,uj). 
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When q is small, we can write the total polarization as 

2 a/3 

P(q,u>) = P^u) + P 2 (q,uj) = I (15) 

irwyoyx — 0/4J 

where x = (hid / / y qa) 2 . Eq. (4) then reads 1 — f3\f?>/(x — 3/4) = 0, where (3 = 2V q /ir'y a, 
and we obtain the solution, 



f, ^3 /, , 4/3 V3 r SV, 

If we compare the result in Eq. (16) with that of Eq. (14), we see that the coefficient in 
front of V q has been doubled. This is because now we have two subbands, while in deriving 
Eq. (14) we assumed there is only one subband. More generally, if there are N (partially) 
occupied, symmetric subbands with degenerate Fermi wavevectors and the same bandwidth, 
we should have the plasmon energy 



/ 2NV 

hu = qhv F Jl + —^. (17) 
V uTivf 

The factor of 2 in front of N can be traced to spin degeneracy. If the subbands are asymmetric 
and each crosses the Fermi level only once, then iV in Eq. (17) should be replaced by N/2. 

From Eqs. (9) and (16), we find the Luttinger parameter g for an armchair carbon 
nanotube to be 

9 = V^l + ^f- ( 18 ) 
For zigzag carbon nanotubes, k y is in the circumferential direction and can only take on 
quantized values k y = v2n/na^ with v = 0, 1,2, ...,n. For v = 2n/3 (remember that n is a 
multiple of 3 for metallic zigzag tubes), we obtain from Eq. (6) the energy dispersions of 
the subbands crossing the Fermi level, 

E hk = ±2 lo sinC^). (19) 

Figure 2 shows the dispersions of the subbands. For the % — 1 subband, the proper integra- 
tion ranges for Eq. (11) are [0,n/^/3a — q] for the first integrand, and [0,7r/\/3a] for the 
second integrand. For the i = 2 subband, the proper integration ranges are [— ir/\/3a, 0] 
for the first integrand, and [— n/^/3a + q,0] for the second integrand. Unlike the armchair 
tubes where the subbands crossing the Fermi level are symmetric, non-degenerate A bands, 
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for zigzag tubes the subbands crossing the Fermi level are asymmetric, doubly degenerate 
E bands. 11 Taking into account the degeneracy, the total polarization is 

_. , N 16 r , Tiw + b — E^taniy/Sqa/IG) 

P(q,u) = — -= [In -= 

a\/?>TiE w hu + b + E w tan(yZqa/l§) 

hu - b + E w tan(\/?>qa/lQ) 

+ in 7= J) (20) 

hu — b — E^tan^y 3qa / 16) 

where b = 47 sm(v / 3<2 , a/8), and E u = ^b 2 — (hu) 2 . For q — > 0, we have 

d/ \ 87o(ga) 2 

a7rV3[|(/i^) 2 - (70^)^ 

From Eqs. (21) and (4), we find that the Luttinger parameter g for zigzag carbon nanotubes 
is the same as that for armchair carbon nanotubes. 

The zigzag tube has four subbands crossing the Fermi level, but each subband crosses the 
Fermi level only once, unlike the armchair tube case where there are two subbands crossing 
the Fermi level but each subband crosses the Fermi level twice. As a result, for zigzag tubes, 
N in Eq. (17) should be replaced by N/2. 

Chiral carbon nanotubes have four asymmetric subbands crossing the Fermi level, and 
each subband crosses the Fermi level only once. Hence N in Eq. (17) should also be 
replaced by N/2. The result is that all single- walled metallic carbon nanotubes have the 
same Luttinger parameter g. 

We have also investigated the effect of changing the Fermi level by applying a gate voltage. 
For an armchair carbon nanotube, lowering the Fermi level would cause the two subbands 
to have different Fermi wavevectors. The i = 1 subband has a smaller Fermi wavevector 
kpi = (2/a)arccos(0.5 — E F /2 r y ), while the i — 2 subband has a larger Fermi wavevector 
kp2 = (2/ a) arccos (0.5 + E F /2 r y Q ). Replacing k F by kpi in Eq. (12) gives the polarization of 
the first subband, while replacing k F by k F 2 in the same equation gives the polarization of 
the second subband. We find that the plasmon energy or the Luttinger parameter g depends 
on the Fermi energy E F very weakly. To second order in E F , we obtain 

»''= 1+ mjgm& (22) 

where go is the g value for E F = 0. To first order in E F , the g value is unaffected by the 
lowering of E F from zero. This can be attributed to the fact that the two subbands have 
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nearly linear dispersions near the K points. Data in Ref. 6 suggest that g = 0.185, and 
(3 = 12.2. Using these values, the above relation can also be written as 

g/g = 1 + 0.107(^) 2 . (23) 
7o 

Band structure calculations have shown that the carbon peapod based on the (10,10) 
carbon nanotube has N = 4 symmetric subbands crossing the Fermi level with nearly 
degenerate Fermi wavevectors. 16 From Eq. (17) we expect the g value for a C 60 @(10,10) 
carbon peapod g p to be smaller than that for the (10,10) carbon nanotube g t . If the four 
subbands of the carbon peapod have the same bandwidth, g p and g t would be related by 

l/gl = 2/gl - 1. (24) 

However, since the calculations in Ref. 16 indicate that two of the subbands have narrower 
bandwidth than the others, the above relation would hold only approximately. Still, experi- 
mental verification of this approximate relation should provide a test for the band structure 
calculations and reveal the connection between Luttinger liquids in the carbon nanotube 
and the carbon peapod. 

The work presented here implies that the RPA is suitable for studying plasmons in carbon 
nanotube systems regardless if the system is a Fermi liquid or Luttinger liquid. Recently, a 
RPA theory has been developed for carbon nanotube bundles. 17 In the case of a 2D electron 
gas in a strong magnetic field, the RPA has been shown to be very good for wavevectors 
up to 2 times the inverse magnetic length, beyond which it is necessary to include exchange 
self-energy and ladder diagrams in a generalized RPA (GRPA), and the GRPA works well for 
arbitrarily large wavevectors. 18 ' 19 We expect similar things to happen in a carbon nanotube 
system. In our future work we will compare the RPA and GRPA for carbon nanotube 
systems to determine the range of validity for RPA alone, and the significance of corrections 
due to exchange self-energy and ladder diagrams. 
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Figure Captions 

Figure 1. For an (n,n) armchair carbon nanotube, only two symmetric, non-degenerate 
A bands cross the Fermi level E F = 0. The energy dispersions of these two subbands are 
given by Eq. (7). The subbands cross each other and the Fermi level at K points (Dirac 
points) k y = ±27r/3. If the Fermi level is lowered to the horizontal dashed line by applying 
a gate voltage, the two subbands will have different Fermi wavevectors. 

Figure 2. For a metallic zigzag carbon nanotube, asymmetric, doubly denegerate E bands 
cross the Fermi level. The energy dispersions of the subbands are given by Eq. (19). 
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